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Abstract 

We investigate the quark angular momentum in a model with the nucleon being a quark and a spectator. Both scalar and axial- 
vector spectators are included. We perform the calculations in the light-cone formalism where the parton concept is well defined. 
We calculate the quark helicity and canonical orbital angular momentum. Then we calculate the gravitational form factors which 
are often related to the kinetic angular momentums, and find that even in a no gauge field model we cannot identify the canonical 
angular momentums with half the sum of gravitational form factors. In addition, we examine the model relation between the orbital 
angular momentum and pretzelosity, and find it is violated in the axial-vector case. 

Keywords: proton spin, orbital angular momentum, gravitational form factor, pretzelosity 


1. Introduction 

Hadrons are bound states of the strong interaction which 
is described by the quantum chromodynamics (QCD) in the 
framework of Yang-Mills gauge field theory. One of the central 
problems in particle physics is to determine nucleon structures 
in terms of quark and gluon degrees of freedom. The decom¬ 
position of the proton spin is one of the most active frontiers 
in recent years. Although the total angular momentum of an 
isolated system is well defined, the decomposition to each con¬ 
stituent of a relativistic composite particle, such as the proton, 
is non-trivial and of great interest. 

The observation that only a small fraction iHi (about 30% 
in recent analysis [3l [4J, |5[]) of the proton spin is carried by 
quark spins has puzzled the physics community for more than 
two decades. This result severely deviates from the naive quark 
model where the proton spin is from quark spins. Many possi¬ 
ble ways to understand the “proton spin crisis” have been pro¬ 
posed, such as to attribute the remaining proton spin to the or¬ 
bital angular momentum (OAM) and/or the gluon helicity. Due 
to the Wigner rotation effect 16[] which relates the spinors in 
different frames, the constituent’s spin of a composite particle 
in the rest frame can be decomposed into a spin part and a non¬ 
vanishing OAM in the infinite momentum frame (IMF) or light- 
cone formalism where the parton language is defined [7, [&j2]. 
Therefore, the OAM plays an important role in understanding 
the “proton spin puzzle”, although the gluon helicity also con¬ 
tributes a large fraction [0. However, the decomposition of 
proton spin, especially the definition of OAM, is still under con¬ 
troversy. 

A most intuitive decomposition is to divide the proton 
spin into quark spin, quark orbit, gluon spin and gluon orbit 
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terms [ IT ]: 

S q + L q +S g + L g = -, (1) 

where the quark orbit operator is defined as 

L q = -i{f/y + r xVifr. (2) 


But the L q , as well as S g and L g , is not obviously gauge- 
invariant and thus renders the physical meanings in common 
situations obscure. To solve this problem, an explicitly gauge- 
invariant decomposition is proposed [ 12]: 


iq+L'q+J'g- -, 


( 3 ) 


where each term is obviously gauge-invariant. It shares the 
same definition for the quark spin operator in £[}. but takes a 
different definition for the quark orbit operator as 


Lq = iij/y + r x Di//, (4) 

where D = -V - igA is the covariant derivative. In this decom¬ 
position, the total angular momentum for each parton flavor is 
usually supposed to be related to the sum of two gravitational 
form factors: 

J q l g = \[A q/g (0) + B q/g m, ( 5 ) 


where the two form factors A and B can be measured through 
the deeply virtual Compton scattering (DVCS) process. 

Recently, Chen et al. revived the idea to decompose the 
gauge potential A M into a pure gauge term ,4 ; r /' rc , which plays 
the role on gauge symmetry and only has to do with unphys¬ 
ical degrees of freedom, and a physical term A„ hy , which in¬ 


volves the two physical degrees of freedom Jl3, 1 14ll . With 
this approach, many more decom posi tion versions were pro¬ 
posed U5, 3. 17]. As observed in [18] and discussed in details 
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in ii;* izutl . this kind of split introduces a so-called Stuekel- 
berg symmetry which copies the group of gauge symmetry but 
acts on the fields on a different manner. Thus, the approach 
of Chen et al. can be viewed as a gauge invariant extension 
(GIE) based on a Stuekelberg symmetry fixing procedure. This 
procedure is essentially a choice of the physical term which is 
frame dependent, and therefore may result in different decom¬ 
position versions which actually correspond to different physi¬ 
cal objects El. Nowadays, all the decompositions are usually 
classified into two groups [22. 231. the canonical version and 
the kinetic (or mechanical) version. Due to the GIE procedure, 
they are both in principle measurable without gauge-invariance 
breaking. 

In this letter, instead of focusing on the controversy on which 
version is more physical, we perform a calculation of quark an¬ 
gular momentums in a spectator model with the nucleon to be 
a struck quark and a spectator. We find that even in this no 
gluon model we cannot identify the canonical angular momen¬ 
tum with half the sum of two gravitational form factors. Then, 
we also examine the model relation between the pretzelosity 
and OAM, and find that the relation is violated in the axial- 
vector case. 


2. Light-cone spectator model 


Hadrons are the eigenstates of the light-cone Hamiltonian 
Hlc = 2P + P~ - P\ with invariant mass square as the eigen¬ 
values. Quantized at a fixed light-cone time r = (t + z)/ V2, 
one may have unambiguous definition on the constituents, and 
hence a hadron state can be expanded on a complete basis of 
Fock states as 112411 


| H,P + ,P ± ,S Z ) = 


N 

x ^ (2) ^ k± j^ n i H l" ; x p 

j =i 


(6) 


where N is the number of constituents of the Fock state [ n), 
Xi, k ± i and A, are the light-cone momentum fraction, intrinsic 
transverse momentum and light-cone helicity carried by the z'-th 
constituent respectively. The i/z„/ W is the light-cone wave func¬ 
tion (LCWF) which decribes the probability amplitude to find 
the Fock state |«) in the hadron state | H). 

In the spectator model, the proton is viewed as a struck quark 
and a spectator which contains the remaining part. Then the 
proton state is expressed as 


\P) 



dxd 2 k ± 

2(27r) 3 V*( 1 — jc) 


fqD(x, k ± )\qD\ X, k j_, d>, 


(7) 


separation. One can effectively introduce the quark-spectator- 
proton vertex in the Lagrangian as [25] 

-C, = g s TV>i/r + g^y^ysA^ + h.c., (8) 


where if/, V F, cf> and A t , are the operators of quark, proton, scalar 
diquark and axial-vector diquark fields. Some suitable form 
factors are included in the effective couplings g s/v to describe 
the structures. Then with the Dirac structure, we write down 
the quark-spectator LCWFs as 


i/s S A(x,k_ l) = 1-. ^t» s (.v. k ± ), 

V2F V2 P+ 

. v a t i \ u(k. A.) , u(P, A) v 


,w(x,k_ l) = 


V2F 


<0, A )/75 


V2F+ 


(9) 

<p\x,k ± -), (10) 


where A, A’ and A are the light-cone helicities of the quark, 
axial-vector spectator and proton, and the superscripts S and 
V denote the type of the spectator. The k , p and P are the 
momentums carried by the quark, spectator and proton respec¬ 
tively. For the Dirac spinor u and polarization vector e, we 
adopt the Lepage-Brodsky convention Ehh . The (p(x. k±) is a 
spin-averaged momentum space wave function. In the follow¬ 
ing calculations, we choose the form as (27, 28] 


*s/v, 


(x, kj = 


&s/ 


tvr 


2 \ 


M 2 - 


Mq + k ^ 
1 - x 


( 11 ) 


where m. Mo and M are the masses of the quark, spectator 
and proton, and g s/v is a coupling constant. This choice corre¬ 
sponds to an effective coupling in the Lagrangian as g sr/ (k 2 ) = 
g s/v /(k 1 - m 2 ). It respects the Lorentz invariance (29] and leads 
to the polynomiality property of generalized parton distribution 
(GPD) moments [27], Replacing the quark mass m by a cut¬ 
off parameter, one can get the form induced by a dipole form 
factor HQ. The LCWFs in ([9]) and (flOl) are normalized as 


' dxd 2 k 

16 7T 3 

dxd 2 k ± 

167T 3 


f 


( 12 ) 


A 2 /v |^ s/v (^,*J| 2 - 1, 


where Ns/v is the normalization factor for spin states: 
, (m + xM) 2 + k 2 , 

N i -3—- 


x■ 

2(1 + x 2 )k\ + 2(1 - x) 2 (m + xM) 2 
x 2 (l - x) 2 

(k\ - xM 2 d - (1 - xfmM ) 2 
x 2 ( 1 - x) 2 M 2 d 
( m + M) 2 k\ 




(13) 


(14) 


x 2 M 2 d 


where x and k L are the light-cone momentum fraction and in¬ 
trinsic transverse momentum carried by the quark. The D rep¬ 
resents the spectator. Constrained by the quantum numbers of 
the quark and proton, the spectator can only be either a scalar or 
an axial-vector, and the axial-vector one is necessary for flavor 


3. Quark angular momentum and gravitational form fac¬ 
tors 

In the quark-spectator state expansion 0, the Fock states 
| qD) are the eigenstates of quark spin S q and total OAM L. 
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Thus, we can express the expected value of quark spin in a po¬ 
larized proton as 


|f s {(x,k ± )| 2 - |^ s }(x,fc L )f 


-u 

1 C* dxd k\ Q 9 Q O 

= 2 J ~[^rWl(x,k ± )N 2 s \4> s (x,k ± )\ 2 


dxd 2 k ± 

167T 3 

dxd 2 k 


(15) 


where 


o (m + xM) 2 - k 2 

W*(x, k ± ) = J ± 


in the scalar case, and 

' dxd 2 k ± 
16tt 3 

dxd 2 k ± 
I6n 2 


C v _ 

~ 


u 

-\I 


(m + xM) 2 + k\ 

^ [|iA V |i'(^ *l)| 2 - \>p v \v(x, k ± )| 2 
wJ(x,k ± )N^\<p v (x,k ± )\ 2 , 


(16) 


(17) 


where 


v 2(1 + x~)k- - 2(1 — x) 2 (m + xM)~ 

Wj (x, k ± ) =- - 


Nlx 2 (\ - x) 2 


(k 


■ xMl - (1 - xfmM) 2 


Nr,x 2 ( 1 - x) 2 M 2 


(18) 


(m + M) 2 k]_ 

N 2 x 2 Mj ) ’ 

in the axial-vector case. The W factors reflect relativistic ef¬ 
fects, such as the Wigner rotation effect. With the LCWFs 
in (fTTb . ijTT} and ( 1 1 71 i are analytically integrable. Substituting 
the quark and spectator mass parameters which are fitted to the 
electromagnetic form factors in [1281 1, we get the numerical val¬ 
ues as 


S* = 0.383, 
Sj = -0.135. 


(19) 

( 20 ) 


Similar to the quark spin, the total OAM can be expressed as 
' dxd 2 k ± 




I6n 2 
dxd 2 k ± 
167r 3 
dxdrk 


W£(x, fc±))Vg|0 s (x, kj.)|", 


-f , 

= f ^^[-\t V l(x,kJI 2 + \^l(x,k ± )l 2 


( 21 ) 


f 


l^ V | 0 fe^)| 2 + 2|^ v |_(x,* ± )| 

dxd 2 l( L 


( 22 ) 


167T 3 

where the W factors are 


wi(x,k ± )N 2 \^(x,k ± )\ z 


The L s and L v are expected values of the total intrinsic OAM 
operator 

d 


L = 


l Yj kl 


i=l 


dk. 


(25) 


Among N intrinsic transverse momentums, only N - 1 of them 
are independent. To get the canonical OAM carried by the 
quark, one may evaluate the expected value of the intrinsic 
OAM operator with respect to the transverse center [3_l[ 321]. 
For the quark-spectator system, the quark intrinsic OAM can 
be evaluated from [33 ] 


i( 1 — x)k l x 


dk 


The numerical values are 


L q = 0.089, 

L v q = 0.017. 

Then the quark total canonical angular momentums are 
Jq = Sq + Lq = 0.472, 


Jq =Sq+Lq =-0.118. 


(26) 


(27) 

(28) 

(29) 

(30) 


Including the angular momentum carried by the spectator, the 
angular momentum sum rule is satisfied in both the scalar and 
the axial-vector cases: 


fS_oS rS rS __j_ 

J ~ ^ q + ^q + ~ 2 ’ 

7 v = s;+i,; + sx + i,s = 1 


(31) 

(32) 


In this no gluon model, these values are expected to be equal to 
the kinetic angular momentums, since the difference between 
the two definitions is an interaction term with the gluon. 

However, with the relation <j5j, it is suggested to obtain the ki¬ 
netic angular momentums for each constituent from the sum of 
two gravitational form factors A(Q 2 ) and B(Q 2 ), which can be 
measured through the DVCS process. Similar to the Dirac and 
Pauli form factors, these two gravitational form factors can be 
calculated from the helicity-conserved and helicity-flip matrix 
elements of the energy momentum tensor current as 134 35] 


<p + q, T ir ++ (0)|P, t> = 20 P + ) 2 A(Q 2 ), 


(33) 


(P + q, T ir ++ (0)|p, I) = 2(P + ) 2 (q m ,q ) B( Q 2 ^ (34 > 

where q 2 = -Q 2 is transfered momentum square. Using the 
Noether theorem j36[] and the effective Lagrangian in ([8]), we 
can derive the canonical energy momentum tensor, which dif¬ 
fers from the Belinfante improved energy momentum tensor by 
an antisymmetric term, as 


Wf(x,k _,.) = 


v (m + MYk , 

Wj(x, k ± ) = , 0 / 

L A 2 x 2 M 2 


(m + xM) 2 + k 2 ± ’ 

2t2 2(1- x 2 )k 2 


Nr,x 2 ( 1 - x) 2 


(23) 

K' =-[<A/a> - - ru 

(35) 

(24) 

T m d v =cPcpd v t - F^d v A p - g^U s + X v ), 

(36) 

rf = - gT-Ci. 

(37) 
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where the subscripts q, D and I denote the quark, spectator 
and interaction parts respectively, and -Cq/s/v represents the La- 
grangian of free quark, scalar and vector fields. Taking the plus- 
plus component of the energy momentum tensor, one can eas¬ 
ily find that the contribution from the interaction term vanishes. 
Thus the contributions from the quark and the spectator are sep¬ 
arated as 

Tq + =^[<Ay + <9> - (<9 + «A)y>], (38) 

T+ + =d + (pd + cf> - F +P d + A p , (39) 


where the expression for the quark park is the same as that in 
QCD with the light-cone gauge A+ = 0. 

After some algebra, the quark parts of these two form factors 
are expressed in terms of the overlap of LCWFs as 


A S q(Q 2 ) 


B l(Q 2 ) 


/ i dxd~k ± ct. at, , , 

l6n3 2 j x ^ a (x ’ k ^\(x, k ± ), 

2 M r 

q 1 - iq 2 J 


(40) 


dxdrk ± 

16 A 


Y xfi S \\x, kj, 

A 

(41) 


for the scalar case, and 
dxdrk 


A V q(Q Z ) 


bAQ 2 ) 


= f ~^~r Y X ^ ] *A' (x, (x, k ± ), (42) 

A,A' 

2 m r 

f - iq 2 J 


dxd 2 k ± 

167T 3 


Y ^aAx, k'jVuAx, k ± ), 

A,A' 

(43) 


for the axial-vector case. The k' ± = k ± + (1 - x)q ± is quark 
intrinsic transverse momentum in the final state. At Q 2 = 0, 
they can be expressed as 


A“(0) = 


*,( 0) = 


I 

f 


dxdrk 


J-u/S 


16;r 3 

dxdrk 


W\(x, k ± )N$x-\<fr(x, A:j_)[ , 


x ii/S 


16tt 3 


Wjjix, k ± )N$x\<fr(x, A:j_)| , 


where 


W*(x, k ± )= 1, 

TJ/S , , , 2M(\ - x)(m + xM) 

Wp(x, k ± ) =- 


(m + xM) 2 + k\ 


and 


Aj(0) - 


<( 0 ) = 


dxd 2 k 


167T 3 


Wg(x, k ± )N y x\A(-x, k ± )\ z 


l6n 3 W A (x,k ± )N y x\<p (x, &_l)| 
dxd 2 k 


where 


(44) 

(45) 

(46) 

(47) 

(48) 

(49) 


K =1 

v 2 M(m + M)(k]_ - xM 2 d - (1 - x) 2 mM ) 
Wb = N 2 x 2 M 2 d 

4 M(m + xM) 

N%~x ' 


(50) 

(51) 


The numerical values are 

A®(0) = 0.290, (52) 

B*(0) = 0.422, (53) 

Aj(0) = 0.294, (54) 

B v q ( 0) = -0.370. (55) 

Including the contributions from the spectator, the momentum 
sum rule: 

A s/V (0) - A^ /v (0) + A® /v (0) - 1, (56) 

and the anomalous gravitomagnetic moment sum rule & 

£ s/v (0) = B^( 0) + B s ^( 0) = 0, (57) 


are satisfied in both cases. 

Therefore, the total angular momentum is equal to half the 
sum of these two form factors 

7 S/V = I[A s/v (0) + B s/V (0)]. (58) 


But, comparing the quark part with the quark canonical angular 
momentums calculated above, we find that even in such a no 
gluon model, where the difference between the canonical and 
kinetic operators L q and L' q has no contributions, one cannot 
identify the canonical angular momentums with half the sum of 
two gravitational form factors for each constituent: 


i[A®(0) + fl*(0)]= 0.356 */ ? s , (59) 

^[Aj(0) + Bj(0)] = -0.038 4 7 9 v . (60) 


In other words, the relation (0) is violated in this model, in both 
the scalar and the axial-vector cases. 

Here we specified a form, the Hwang-Mueller prescrip¬ 
tion [127, 28], for the spin-averaged LCWFs (p s ^ v (x,k ± ) to get 
the numerical values quantitatively. This form respects the 
Lorentz invariance and produces the GPD polynomiality prop¬ 
erty. Technically, it makes all the integrals in this letter an¬ 
alytically integrable. Apart from this prescription, there are 
many other choices, such as the Brodsky-Huang-Lepage (BHL) 
prescription M 39 ; 0, the Terentev-Karmanov (TK) pre¬ 
scription 0, 14211 . the Chung-Coester-Polyzou (CCP) prescrip¬ 
tion |43|. the Vega-Schmidt-Gutsche-Lyubovitskij (VSGL) pre¬ 
scription 0 and so on. These prescriptions have nothing to 
do with any W factors which are determined by Dirac struc¬ 
tures. Therefore, any choice will not change any conclusion in 
this letter, although it will indeed change the numerical values 
quantitatively. 


4. Pretzelosity and orbital angular momentum 

The pretzelosity, denoted as hf T , is one of the eight lead¬ 
ing twist transverse momentum dependent parton distributions 
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(TMDs). It represents the probability to find a transverse po¬ 
larized quark in a perpendicularly transverse polarized pro¬ 
ton, and can be measured through the single spin asymme¬ 
tries in the semi-inclusive deep inelastic scattering (SIDIS) pro¬ 
cess |45.46]. Based on some model calculations, it is suggested 
to relate the pretzelosity to the OAM 1471, 32]. 

In the scalar case, the pretzelosity is expressed as 

2 M 2 


hfr(x,k ± ) = 


1 

167r 3 


-\(f> (x, &j_)|“ 


where 


N 2 2M 2 ( 1 - A') 3 
167t 3 [k\ + A 2 ] 4 


hr = xM~ d + (1 - x )m~ — jc( l — x)M~ 


(61) 


(62) 


Its first transverse momentum moment is equal to the OAM of 
the quark-spectator system with a minus sign 0: 


L s = 


f 


dxd 2 k i 


2 M 2 


4 hfS(x, k ± ) = 0.117, 


(63) 


which is also equal to the difference between the helicity and 
transversity |9j,|25(]. Similarly, the intrinsic OAM carried by the 
quark can be expressed with the pretzelosity as 




-/ 


dxd 2 k A 


(1 - x)hfj(x,k ± ). 


2 M 2 


(64) 


However, in the axial-vector case, the pretzelosity is ex¬ 
pressed as 

2 M 2 (m + M) 2 


hfj (x, k ± ) = 


1 

167T 3 

N 2 


I <f> V (x, k _ L )| 2 
2M 2 (m + M) 2 (l - x) 3 


x 2 M 2 d 


(65) 


167r 3 + A 2 ] 4 

Its first transverse momentum moment is 

' dxd 2 k ± 


f dxd2k ^^ h tr = - f ^|^ V loOaj [ 2 


( 66 ) 


which is equal to the third term of the L v in (l22l > with a minus 
sign. Hence, simple extensions of the relations (l63l > and ( 164b to 
the axial-vector case are not justified, because in this case the 
model is not spherically symmetric as demonstrated in [32, 48]. 
If only the transverse polarizations for the axial-vector are in¬ 
cluded, the pretzelosity will vanish in this model, but the in¬ 
trinsic OAM is still non-vanishing. Therefore, without further 
assumptions, there is no general relations between the OAM 
and pretzelosity. This is consistent with our intuition, since the 
OAM is essentially a correlation between the coordinate and 
the momentum, while the pretzelosity, as a leading twist TMD, 
only contains the momentum information. If the term in (l66l) 
plays the dominant role in L v , the relations ( I6.H and ( 164b might 
be extended to the axial-vector case as approximate relations 
with some correction factors phenomenologically. Comparing 
([65l> with the OAM expression in (l22l >. this correction factor is 
written as 


C v (x, k ± ) = 


(1 - x)(m + MY 


■2(1 +x )M 2 d 


(1 — x)(m + M) 2 


(67) 


5. Conclusions 


In this letter, we investigate the quark angular momentum in 
a spectator model. The calculations are performed in the light- 
cone formalism where the parton concept is well defined. 


Nowadays, there are many decomposition versions for the 
proton spin. All these decomposition versions are usually clas¬ 
sified into two groups, the canonical version and the kinetic ver¬ 
sion, and the main difference between them is the definition of 
the OAM [22[ 23]. With the GIE procedure, both of them are 
in principle measurable without gauge-invariance breaking. In 
this study, regardless of the dispute on which one is more phys¬ 
ical, we perform the calculations in a no gluon model where the 
results from these two definitions are expected to be the same. 


Considering the Dirac structure, we write down the LCWFs 
of the quark-spectator system with the scalar and axial-vector 
couplings. Then we calculate the spin and intrinsic canonical 
OAM carried by the quark in both cases. By including the con¬ 
tributions from the spectator, the angular momentum sum rule is 
satisfied. Taking the relation ([5} as an assumption fl3l . we cal¬ 
culate the so-called kinetic quark angular momentum through 
two gravitational form factors. As a direct result of the momen¬ 
tum fraction and the anomalous gravitomagnetic moment sum 
rules, the total angular momentum is equal to half the sum of 
these two form factors. However, for each constituent, we can¬ 
not identify half the sum of the form factors with the canonical 
angular momentums, which are expected to be the same as the 
kinetic ones in a no gluon model, in either the scalar case or 
the axial-vector case. In other words, the relation (0) is vio¬ 
lated in this model, even though no gluon degrees of freedom 
are introduced. 


In principle, the proton spin decomposition should be under¬ 
stood with explicit calculations in QCD. But due to the nonper- 
turbative nature of QCD at hadron scale, it is almost impossible 
to perform such an example at present. Instead, an example in 
QED, which has very similar operator structure to that in QCD 
up to a color factor, was performed recently, but surprisingly 
the relation (0 failed to match order by order in perturbative 
calculations J49]. An explicit calculation in this letter indicates 
that neither does the “general” relation between the kinetic an¬ 
gular momentum and the gravitational form factors satisfy in 
the spectator model, which has been widely used to investigate 
the nucleon structure phenomenologically. Therefore we need 
more careful scrutiny concerning issues related to the orbital 
angular momentum of a composite system. 
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Appendix A. Axial-vector with only transverse polariza¬ 
tions 


If only transverse polarizations, i.e. A' = +1, of the axial- 
vector spectator are included, the spin states normalization fac¬ 
tor in (IT4l) is expressed as 


Ny = 


2(1 + x 2 )k 2 ± + 2(1 - x) 2 (m + xM) 2 
x 2 (l - x) 2 


(A.l) 


Such situation applies when the spectator to be massless. Then 
the W factors are correspondingly written as 


W^(x,k ± ) 

Wj(x,k ± ) 

Wj(x,k ± ) 

Wft(x, k ± ) 


(1 + x 2 )k]_ - (1 - x) 2 (m + xM) 2 
(1 + x 2 )k]_ + (1 - x) 2 {m + xM) 2 ’ 

(1 - x 2 )k 2 ± 

(1 + x 2 )k 2 ± + (1 - x) 2 (m + xM) 2 

= 1 , 

2 M(m + xM)x( 1 - x) 2 
(1 + x 2 )k 2 ± + (1 - x) 2 {m + xM) 2 


(A.2) 

(A.3) 
(A.4) 
(A.5) 
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